ON 

On 
On 



Boundary Kondo impurities in the generalized 
supersymmetric t — J model 

Heng Fan"' 6 , Miki Wadati a , Rui-hong Yue 6 

qj . a Department of Physics, Graduate School of Science, 

j_i | University of Tokyo, Hongo 7-3-1, 

c/2 \ Bunkyo-ku, Tokyo 113-0033, Japan. 

^Institute of Modern Physics, P.O.Box 105, 
Northwest University, Xi'an 710069, P.R.China. 



Ctf 



-O '. February 1, 2008 

C 
O 

o 



(N 



Abstract 

We study the generalized supersymmetric t — J model with Kondo impurities in the boundaries. 



We first construct the higher spin operator K-matrix for the XXZ Heisenberg chain. Setting the 
g\ ■ boundary parameter to be a special value, we find a higher spin reflecting K-matrix for the super- 

symmetric t — J model. By using the Quantum Inverse Scattering Method, we obtain the eigenvalue 
and the corresponding Bethe ansatz equations. 
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1 Introduction 

There has been extensive interests in the investigation of low-dimensional correlated electron systems with 
impurities. Recently, using renormlization group techniques, Kane and Fisher p) studied the transport 
properties of a ID interacting electron gas in the presence of a potential barrier. They showed that a single 
potential scatter may dramatically influence the physics in the presence of repulsive e — e interactions. The 
system behaves like a Tomonaga-Luttinger liquid rather that a Fermi liquid. Some different techniques 
were also applied to study similar systems B H]. The Kondo impurities in a Tomonaga-Luttinger liquid 
have been investigated in great detail ||, ||, ||, [TJ] . 

Attempts to study the effects due to the presence of impurities in ID quantum chains in the framework 
of integrable models have a long succesful history|q|-|[l4|. Andrei and Johannessonj| studied an arbitrary 
spin S embeded in a spin-f/2 Heisenberg chain. This method was generalized to other cases. Recently 
the supersymmetric t — J model with impurities has attracted considerable interests. The Hamiltonian 
of the t — J model includes the near-neighbour hopping (t) and antiferromagnetic exchange (J) |L5L |l6| 



H = Z^ 



j = l I <T=±\ ' ) 



It is known that this model is supersymmetric and integrable for J = ±2t |Ia Hq]. The supersymmetric 
t — J model was also studied in Refs.]!^, |20| , pl| , g2|, for a review, see Ref.(g3| and the references therein. 
Essler and Korepin et al [p2| showed that the one-dimensional Hamiltonian can be obtained from the 
transfer matrix of the two-dimensional supersymmetric exactly solvable lattice model |2l], |2J] . 

By use of the Qunatum Inverse Scattering Method (QISM) fl25|, the supersymmetric t — J model 
with higher spin impurity was first investigated in the periodic boundary conditions p4| . Recently, the 
supersymmetric t — J model with impurities have been studied extensively in both periodic and reflecting 
(open) boundary conditions [g6|, g7|, |28|, p^| . 

The open boundary condition was studied extensivly in the last decade. There have been several 
methods to study the the problem of open boundary condition |3fj, £51j. In the end of 80's, Sklyaninf32| 
proposed a systematic approach to handle the open boundary condition problem in the framework of 
the QISM. Besides the Yang-Baxter equationpjl, the reflection equation proposed by CherednikBJ] also 
plays a key role in proving the commutativity of the trasfer matrix. We know that the Hamiltonian of 
the model is usually written as the logarithmic derivative of a transfer matrix at zero spectral parameter. 
The boundary terms in the Hamiltonina are determined by the reflecing K-matrix which is a solution 
to the reflection equation. In the usual boundary problem, the K-matrix is a c-number matrix. The 
operator K-matrices which determine the Kondo impurities in the Hamiltonian are studied recently for 
several models |B5), including the supersymmetric t — J model [E6] E7J. 

The Hamiltonian (1) of the supersymmetric t — J model can be obtained from the transfer matrix 
constructed by the rational R-matrix. We can also use the trigonometric R-matrix to formulate the 
transfer matrix. The corresponding Hamiltonina is a generalization of the original supersymmetric t — J 
model[Q. This Hamiltonian satisfies a symmetry of the quantum group SU q (2\l). In this paper, we shall 
study the generalized supersymmetric t — J model with higher spin boundary impurities. The operator 
K-matrix is first constrcuted for the XXZ Heisenberg spin chain with higher spin impurities. We then 
find a higher spin operator K-matrix for the supersymmetric t — J model. Using the graded algebraic 



Bcthe ansatz method, We obtain the eigenvalue of the transfer matrix and the Bethe ansatz equations. 
The paper is organized as follows: We introduce the model in section 2. In section 3, we study the 
XXZ spin chain with higher spin Kondo impurities and present the higher spin reflecting matrices for 
the generalized supersymmetric t — J model are obtained. In section 4, using the nested algebraic Bcthe 
ansatz method, we obtain the eigenvalues of the transfer matrix for the generalized supersymmetric t—J 
model. Section 5 includes a brief summary and discussions. 

2 The Model 

We first review the generalized supersymmetric t — J model. For convenience, we choose similar notations 
as those in |22| and our previous paper [B7l. The Hamiltonian of the generalized supersymmetric t—J 
model take the following form: 

N 

H = J2Y1 HA 1 ~ n J;-<^) c 3+iA 1 - nj+i,-<r) + c ]+iA l - nj+i,-a)c 3 +iA l - nj,-a)] 



- 2 D!(^+i + s .Ai) + «"fa)^s;+i - - s(n} 



N 



+isin( V ) 2JS!n j+1 - S z 0+1 n 3 ]. (2) 

i=i 
When the anisotropic parameter rj = 0, this Hamiltonin reduces to an equivalent form of the original 
Hamiltonian (1). The operators Cj >a and c- a mean the annihilation and creation operators of electron 
with spin a on a lattice site j, and we assume the total number of lattice sites is N, a = ± represent spin 
down and up, respectively. These operators are canonical Fermi operators satisfying anticommutation 
relations 

{4,<7' C .7>t} = <M<7T- (3) 

We denote by rij_ a — c'- a Cj_ a the number operator for the electron on a site j with spin a, and by 
Uj = J2a=± n h° the number operator for the electron on a site j. The Fock vacuum state |0 > is defined 
as Cj, CT |0 >= 0. Due to the exclusion of double occupancy, there are altogether three possible electronic 
states at a given lattice site j 

|0>, |T> J -c], 1 |0>, \i> j =c] < _ 1 \Q>. (4) 

Sj, Sj, Sj are spin operators satisfying su(2) algebra and can be expressed as: 

Sj = 4iC i: _i, 5J = c] _ lCj , , S? = -(n jA - rij-i). (5) 

The above Hamiltonian can be obtained from the logarithmic derivative of the transfer matrix at zero 
spectral parameter. In the framework of QISM, the transfer matrix is constructed by the trigonometric 
R- matrix of the Perk-Schultz model |3q| . The non-zero entries of the R-matrix are given by 

R(X) a a a a = sin(v + e a X), 

RWtb = (-ir^szn(X), 

R(\)f a = e ts "^ a -V x sin( V ),a^b, (6) 



where e a is the Grassman parity, e a = for boson and e a = 1 for fermion, and 

/ j\ f 1: if a > 6 
s*gn(a-b) = ! [ ^ . f q < ^ 

This R-matrix of the Perk-Schultz model satisfies the usual Yang-Baxter equation: 
7? 12 (A - [*)Ri 3 (\)RM = R23(n)Ri 3 WRi2(X - n) 



(7) 



(8) 



In this paper, we shall concentrate our discussion only to the Fermionic, Fermionic and Bosonic case 
(FFB), that means ei = e 2 = 1, £3 = 0. And we shall use the graded formulae to study this model. For 
supersymmetric t — J model, the spin of the electrons and the charge 'hole' degrees of freedom play a 
very similar role forming a graded superalgebra with two fermions and one boson. The holes obey boson 
commutation relations, while the spinons are fermions p3. The graded approach has an advantage of 
making clear distinction between bosonic and fermionic degrees of freedom |39[| . 

Introducing a diagonal matrix H a d c = (— ) Ca£c <5ah^cd, we change the original R-matrix to the following 
form 



R(X) = ILR(A). 

From the non-zero elements of the R-matrix R°£, we see that e a + £& + e c 
the R-matrix satisfies the graded Yang-Baxter equation 



(9) 
Ed = 0. One can show that 

(10) 



Explicitly the R-matrix is written as 



R(X) = 



( o(A) 
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(11) 



where 



a(A) = sin(X — 77), w(\) — sin(X + 77), 6(A) = sin(X), c±(X) — e sin(j]) 



In the framework of the QISM, we can construct the L operator from the R-matrix as: 

L n (X) = 



6(A) - (6(A) - o(A)K x -c_(A)eJ x c_(A)e^ 

- c+ (A)e 1 1 2 6(A)- (6(A) -a(A))e^ 2 c_(A)e^ 2 

c+ (X)e n 13 c + (A)e^ 3 6(A) - (6(A) - w(X))e% 3 



Here e™ & acts on the n-th quantum space. Thus we have the (graded) Yang-Baxter relation 

i? 12 (A - M )£i(A)i 2 (/i) = i 2 (/j)ii(A)i? 12 (A - M ). 



(12) 



(13) 



(14) 



Here the tensor product is in the sense of super tensor product defined as 

(F®G) b a d c = F b a G d c (-)^+^. (15) 

Except in Section 3.1, all tensor products in this paper are in the super sense. 

The row-to-row monodromy matrix T/v(A) is defined as a matrix product over the N operators on all 
sites of the lattice, 

T„(A) - L aN (X)L aN _ 1 (X) ■ ■ ■ L al (X), (16) 

where the subscript a represents the auxiliary space, and the tensor product is in the graded sense. 
Explicitely we write 



{[TW] } 



1 b \ 

01---0N 

= iiv(A)^^_i(A)^-^----£ 1 (A)^ 1 (-l)Sr= 2 K-+^)E:: i 1 ^ (17) 

By repeatedly using the Yang-Baxter relation (|l4[), one can prove easily that the monodromy matrix also 
satisfies the Yang-Baxter relation 

R(X - ^)Ti(A)T a (/i) = T a Gu)Ti(A)fl(A - //). (18) 

For periodic boundary condition, the transfer matrix T per i(X) of this model is defined as the supertrace 
of the monodromy matrix in the auxiliary space 

r P eri(A) = strT(X) = ^(-l) e "T(A) aa . (19) 

As a consequence of the Yang-Baxter relation (nq) and the unitarity property of the R-matrix, we can 
prove that the transfer matrix commutes with each other for different spectral parameters, 

[Vpert (A), T pe ri (/■«)] = 0. (20) 

In this sense we say that the model is integrable. Expanding the transfer matrix in the powers of A, we 
can find conserved quantitcs. The first non-trivial conserved equantity is the Hamiltonian (1), 

H = sin{ji) dH ^ X)] \ x=0 = J2 H Jd+i = E^+i^+iCO), ( 21 ) 



■>bd 



where Pij is the graded permutation operator expressed as P™ = 5 a dSbc(— l) c ° £c . 

In this paper, we consider the reflecting boundary condition case. In addition to the Yang-Baxter 
equation, a reflection equation should be used in proving the commutativity of the transfer matrix with 
boundaries. The reflection equation takes the form p3 

R 12 (X - n)Ki(X)Rn{X + y)K 2 {^) = K 2 {^)R 12 {X + l i)K 1 {X)R 21 (X - fi). (22) 

For the graded case, the reflection equation remains the same as the above form. We only need to change 
the usual tensor product to the graded tensor product. We write it explicitly as 

J2(A-/i)^ a A'(A)giJ(A + /i)g*A'( At )*(-)^+ e 'x)^ 

= ^(/i)^i?(A + M )^ 2 K(A)^i?(A- M )^(-)K+^K 2 . (23 ) 



Instead of the monodromy matrix T(X) for periodic boundary conditions, we consider the double-row 
monodromy matrix 

T(X) = T(\)K(\)T-\-\) (24) 

for the reflecting boundary conditions. Using the Yang-Baxter relation, and considering the boundary 
K-matrix which satisfies the reflection equation, one can prove that the double-row monodromy matrix 
T(A) also satisfies the reflection equation 

R(X- ^l\tT(X)tlR(X + ^)tlfT(^(-)^+^^ 
= T( Ai )^i?(A + / i)^T(A)-i?(A- / i)^(-) ( ^ +e - )£ -. (25) 

Next, we study the properties of the R-matrix. We define the super-transposition st as 

(A s % = A rt (-1)^+ 1 ^ . (26) 

For FFB grading used in this paper, t\ = e 2 = 1, £3 = 0, we can rewrite the above relation explicitly as 

An A 12 B 1 \ st ( A n A 21 C x 

A21 A 22 B 2 = A l2 A 22 C 2 I . (27) 

C x C 2 D ) \ -B ± -B 2 D 

We also define the inverse of the super-transposition st as {A st } st = A. 



One can prove directly that the R-matrix (11) satisfy the following unitarity and cross-unitarity 
relations: 

Ri 2 (X)R 21 (-X) = p(X)-id. 1 p(X) = sin(r) + X)sin(r]-X), (28) 

R?2 (V ~ X)M 1 R 21 1 (A)Mf l = /5(A) • id., p{X) = sin{X)sin{r, - A). (29) 

Here the matrix M — diag.(e 2lTI , 1, 1) is determined by the R-matrix. The cross-unitarity relation can 
also be written as the following form 

{M 1 1 Rir t2 (r 1 -X)M 1 f 2 R^(X) = p(X), (30) 

R^{X){M 1 Rf 1 st H^~X)M l 1 f 2 = p(X). (31) 

In order to construct the commuting transfer matrix with boundaries, besides the reflection equation, we 
need the dual reflection equation. In general, the dual reflection equation which depends on the unitarity 
and cross- unitrarity relations of the R-matrix takes different forms for different models. For the models 
considered in this paper, we can write the dual reflection equation in the following form: 

R^ (p - X)K+ Stl {X)M 1 l Rf 2 ^ {rj-X- p)M 1 K+ st2 (p) 
= K+ st2 (p)M 1 R s 2 t 1 st ^(r 1 - X - p)M 1 - 1 K+{X)R^ st2 {p ~ A). (32) 

Then the transfer matrix with boundaries is defined as: 

t(X) = strK+(X)T{X). (33) 



The commutativity of t(X) can be proved by using unitarity and cross-unitarity relations, reflection 
equation and the dual reflection equation. The detailed proof of the commuting transfer matrix with 
boundaries for super (graded) case can be found, for instance, in Ref.|40| [ll| [l2| etc.. With a normalization 
K (0) = id., the Hamiltonian can be obtained as 



1 rflnf(A) , 
H = -sm(T)) - | A =o 



2 w/ dX 

N-l 



£, p^l^m + Un MKm + s ' r " A :T+;> (0) - <*> 



sir a if+(0) 

3 Higher spin solution to the reflection equation for supersym- 
metric t — J model 

In order to find the higher spin solution to the reflection equation for the generalized supersymmetric 
t — J model, we first construct the higher spin reflecting matrix for the XXZ Heisenberg chain. 

3.1 XXZ Heisenberg chain with higher spin boundary impurities 

The higher spin R-matrix can be constructed by using the fusion procedure E3 . The Hamiltonian of the 
XXZ Heisenberg chain is written as 

N 1 

i=i 

Here a ± = l/2(a x ±a v ) and a x ,a y and a z are Pauli matrices. The R-matrix is known to be the standard 
six-vertex model, 



ria(A) 



/ sin(\ + r)) 

sin(X) sin{rf) 

sin{rf) sin(\) 

\ sin(X + r)) ) 



(36) 



In the framework of QISM, the L-operator constructed by the r-matrix is written as: 

T ,m- ( sin{X+±ri + ±r)a z k ) sin(r))a k \ . . 

Lak[A) - { sin( V )a+ sm(X + \r, - \ W %) J' [6 ' } 

where a represents auxiliary space. As usual, we can construct the row-to-row monodromy matrix 
T a (X) = LaN(X) ■ ■ ■ L a i (A), and we have the Yang-Baxter relation 

r 12 (A - n)Ti(X)T a (jj,) = T 2 ( M )T 1 (A)r 12 (A - M ), (38) 

where the tensor product is a non-graded one. 

Next, we consider the higher spin operators. Let the higher spin L operator take the form J43| |44j 

, { sin{X + S z r]) sin(r))S- \ , . 

[ ' ~ I sin(r])S+ sin(X - S z r/) J ' [ ' 



where S z , S and S' are spin-s operators satisfying the following commutation relations, 

[S Z ,S ± ]=±S ± , [S+,S-] = Sm(2SZ x 7?) . (40) 

sin(rj) 

We also have the following relations for spin-s operator: 

sin(S z 7])sin(rj + S z n) + sin (?7)S~S + 
= sin 2 (r))S + S~ + sin(S z ?])sin(S z n — n) — sin(sn)sin(srj + n) (41) 

A more general relations can be written as 

sin(X + S z ?y)sm(7y + S z ?7 — A) + sin (?7)S~S + 
= sin (?7)S + S~ + sin(X — S z i])sin( — A — S z ?; + n) = sin{\ + sn)sin(sr] + 77 — A). (42) 

One can prove that the higher spin L operator also satisfy the Yang-Baxter relation 

r 12 (A - m)A(A)£ 2 (m) = £a(M)A(A)ri 2 (A - m)- (43) 

Now, Let us consider the reflecting boundary condition. We can find a c-number solution to the 
reflection equation K r (X) — diag. (sin(£ + \),sin(£ — A)), where £ is an arbitrary parameter. This is 
a general c-number diagonal solution to the reflection equation. In particular, if £ — ► — ioo, we find 
K{\) = diag.(e 2lX , 1) is a solution to the reflection equation. 

It is interesting to find a higher spin operator K-matrix. We can construct the operator K-matrix by 
fcxxz(X) = £(A + c)K c (\)£~ 1 (—\ + c), one can find easily that /C(A) is an operator reflecting matrix 
satisfying the reflection equation. Explicilty, the higher spin reflecting K, has the form /Cxxz(A) = 

with 



K(X)\ K{\)\ 
K{\)\ K{\)1 

K-Wi = sin(X — £)sin(\ + c + sn)sin(X + c — n — sn) 

+sin(2X)sin(X + c + S z ri)sin(t; — c + n + S 2 ??), 

K-Wz = —sin(t; + X)sin(X + c+sri)sin(X + c — n — sri) 
+sin(2X)sin(X + c — S z i])sin(^ + c — n + S z ?7), 

K{X)\ = sin(r])sin(2X)sin(^ + c+S z j 1 )S-, 

K{X)\ = sin(r])sin(2X)sin(C~c+S z r 1 )S + . (44) 

By use of the cross-unitarity relation of the r-matrix, the operator reflecting matrix to the dual reflection 
equation can also be found. The eigenvalues of the transfer matrix can be obtained by applying the 
algebraic Bethe ansatz method. These results will be presented in a seperate paper |4q| . 

3.2 Higher spin reflecting matrix for the supersymmetric t — J model 

We know that the generalized supersymmetric t — J model has a SU q (2) symmetry. We suppose that the 
operator K-matrix takes the following form: 

A(X) B(X) 
K(X) = J B(X) C(X) J (45) 

1 



Inserting this matrix into the reflection equation (E3), we can find the following non-trivial relations: 

f (A - ^XK{\)l\f{\ + tillfK^tl = K(») b a lr(\ + »tZK{\)i\r(\ - „)** , (46) 



and 



K(\% 1 K<jj)%=KOj)* l K(\)il, 



(47) 



S aidl sin(X - M)e" i(A+/j) + sin(X + ^ x -^K{\) d a \ 
= e-^-^sin(X + n)K(n)i\ + e^ x +^K^) b a \K(X) d b l, 



(48) 



where all indices take values 1,2, and we have introduced 



ri 2 (A) = 



/ sin(X-ri) 

sin(X) — sin(ri)e~ lX 

—sin(i])e tX sin(X) 

\ sin(X-r)) ) 



(49) 



This matrix f (A) can be obtained from the matrix ( p6|) by a gauge transformation and with a change 
Tj -> -77. Correspondingly, we can show that A(X) = f{X)e^ 2iX K{X)\ 1 B(X) = f(X) e - iX K(X) 2 , C{X) = 
f{X)e- iX K(X)l 1 D(X) = f(X)K{X) 2 satisfy relation @. Substitute these results into relations @||), 
and after some tedious calculations, we find that if we take £ — * —zoo, and /(A) = — 1 / e 2lX sin(X — c — 
r\ — sr\)sin{X — c + sr]), all relations obtained from the reflection equation can be satisfied. So, we finally 
find the higher spin reflecting matrix as 

A(X) = g(X) (e- ilX sin(X + c - sr))sin(X + c + rj + sr]) - sin(2X)sin(u + c - s z 7 ? )e~ l(3A+c+r ' +SSi '' ) ) , 

B(X) = g(X)sin( V )sni(2X)e- li2X - c+s * ,,} S-, 
C(X) = g(X)sin(r 1 )sin(2X)e~ z{2X+c+s * 7 > ) S + , 
D(X) = g(X)(sin(X + c-sri)sin{X + c + T ] + ST 1 )-sin(2X)sin(X + c+S z ri)e-^ x - c - r ' +s ^" ) ), (50) 



where g(X) = l/sin(X — c — rj — sr])sin(X — c + sr]). 

Next, let us consider the higher spin reflecting matrix to the dual reflection equation (B2T). We 
suppose K + has the similar form as K. By direct calculation, we can find Rl 2 lSt2 (X) = I\R 2 i(X)I\ 
with / = diag(—l, —1, 1). For the form ((45[), we have IK(X)I = K(X). Then with the help of property 
[M1M2, R(X)] — 0, we can write the dual reflection equation as 



R12U1 - X)K+ Stl (X)M- 1 R 21 (t 1 - A - ]i)K+ st2 {^M- 1 = 
K+ St2 (^M^R 12 ( V - A - ri)K+ Stl (X)M^R 21 (ri - A). 



We see that there is an isomorphism between K and K + 

K(X) 



K +st (X) = isT(|-A)M. 



(51) 



(52) 



Given a solution to the reflection equation (123) , we can also find a solution to the dual reflection equation 
(pl|). Remark that in the sense of the transfer matrix, the reflection equation and the dual reflection 



equation are independent of each other. We can write the higher spin reflecting matrix to the dual 
reflection equation as 

/ A+(X) B+(X) \ 
K+(X) = B+(\) C+(X) , (53) 

V 1/ 

with 

A + (X) = g + (\)[e 4iX sin(X + c-j] + srj)sin(X + c-2j]-sri). 

-sin(2X - rj)sin{u + £~r/- S z ry)e l(3A+£ -'' +s ~ z '' ) ], 
B+(X) = -g + (X) S in(T 1 )sin(2X-T 1 )e ti2X+£+ ^ + ^ ) S- 1 
C+(X) = -g + {X)sin{7 1 )sin{2X~7 1 )e l{2X - £i -^ + ^ ) S+ 1 
D + (X) = g + (X)[sin(X + c — rj + srj)sin(X + c — 2r/ — st]) 

-sin(2X - r))sin(X + c - rj + S~ z ?7)e i(A - a+ '' +sS7) ], (54) 

where g + (X) — 1/ sin{X — c + r] + sr])sin(X — c — stj). 

Thus we find the higher spin reflecting matrices for the generalized supersymmetric t — J model. We 
should remark that these higher spin reflecting matrices are the kind of 'singular' matrices. It can not 
be constructed directly by the Sklyanin's 'dressing' procedure. In the rational limit, it reduces to the 
result obtained in [p7| . The rational higher spin K-matrix has been analyzied in detail by the projecting 
method |2S|. Our result should also be obtained by the projecting method. 

4 Algebraic Bethe ansatz method for the generalized supersym- 
metric t — J model with higher spin impurities 

4.1 First level algebraic Bethe ansatz 

We denote the double-row monodromy matrix as 

/ .An (A) A 12 (X) Si (A) \ 
T(A) = A 2 i(X) A 22 (X) B 2 (X) . (55) 

V d(A) C 2 (X) V{X) ) 

For later discussions, we introduce the following transformations 

A ab (X) = A ab (X) + 5 ab e 2 ] X * m{7] \ v(X). (56) 

sin(2X + rj) 

As mentioned in section 2, the double-row monodromy matrix satisfies the reflection equation (p5l), we 
have the following commutation relations: 

tWA - u) d2dl 

C dl (X)CM = 12 N 2 C1 . CMCcAV, (57) 

sm{X — fi + rj) 

svn(A + /! + 7])sin(\ — /i) 
S in(2»)sin { ,)e^ ^^ _ ^^ ^ (5g) 



sin(X — /i)sm(2/i -\-rj) ' sin(X + fi + r/) 

10 



fi 2 (A + n + ry)oi 6 c 2 2 r 2 i (A - a*)^, 2 

-4aidi(A)C d2 (/i) = . ,, , , 1 s . 77 r J -^-C C2 (^)A 1 6 1 (A) 

sin{\ + fi + j])sin(\ — fi) 

+ «n(A - p)«n(2A + q) r " (2A + ^^ (A) ^* (m) 

S? n(A + p + V)^n(2\ + q) a m(2 M + q) riz(2A + ^^( A W- ( 59 ) 



Here the indices take values 1,2, and the matrix f is defined in fl49|). 

We define a reference state in the n-th quantum space as |0 >,,= (0, 0. 1)', and reference states for the 
boundary opeators as S"|0 > r = 0,S Z |0 > r = -s|0 > r ,S+|0 > r ^ 0, and S~|0 > ; = 0,S Z |0 >,= -s|0 >/ 
, S+|0 >iy^ 0. The vacuum state is then defined as |0 >= |0 >i <^^ =1 \0 >k <8>|0 > r . Acting the double-row 
monodromy matrix on this vacuum state, we have 

B a {\)\0> = 0, 
C a (A)|0> ? 0, 
D(A)|0> = sin 2N {X + T 1 )\0>, 

A 6 (A)|0> = sm 2N \X)[K{X) b a - 5 ab Sm{ ^^\ ]\Q >= W ab {X)sin 2N (A)|0 >, (60) 



where 



sin{2\ + rj) 



Wi 2 (A) = 0, W 2 i(A) = C(A), 

W n {\) = ff(A) et "^^ 2A \ [e-^ 4A+2 ")am(A + c + e - «j)«n(A + c + 2r? + «,) 

- sm(2A + r])sin{\ + c + r/ + sJ7)e _i(3A+c+3 ' ) - s ' 7) ], 
W fn „-2a sm(2A)sm(A + c + 77 - 577) 

^ 22(A) - C sin(2X + r,) 3 in(X-c + 3 V Y (61) 



The transfer matrix (B3h can be written as 

t(X) = -K + (X) a b A ab (X) + V(X) 

= -K+(X)» a A ba (X) + (l - "$2\ + v) [A+W + D+(X)] ) V{X) - (62) 

Acting this transfer matrix on the ansatz of the eigenvector 

C dl {^i)C d2 {^ 2 ) ■ ■■C dn (ix n )\0 > F dl - d ", (63) 

where F x " n is a function of the spectral parameters /ij, we have 

t(X)C dl ( Pl )C d2 (/i 2 ) • • • C dn (/i„)|0 > F dl - d " 
sin(2X — r/)sin(X — c + r\ — sr])sin(X — c + 2r\ + srj) 
sin(2X + rj)sin(X — c + 77 + srj)sin(X — c — srf) 



xsin 2N 



(A + r,) ft ™£ I f } f f." f " % M ■ -CM10 > F*"4> 
f_ x sin(X + Hi + r))sin{X — /i^) 
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1 

+"»"'<*> n 8i -» (A _ „)» i , i ( A + ,„ + / ■■ '"■» ■ • •c-o-)'"' wi;v.ai° > **"* 

+«.*., (64) 

where u.t. means the unwanted terms, and £W(A) is the so-called nested transfer matrix which can be 
defined, with the help of the relation (B9|) , as 



t {1) WZ- C dl = - K+ Wt {r(A + Mi + r?)£fr(A + /i 2 + vTaltl ■ ■ -r(A + /ii + j/)*^ } 

x^a„ b „(A) {f 21 (A - /*„fc„ ld " • • ■ f2i(A - M2)^ 2 2 r 21 (A - Ml )^ } . (65) 

We find that this nested transfer matrix can be regarded as a transfer matrix with reflecting boundary 
conditions corresponding to the anisotropic case 

^\X)=strK^ + {X')T^{X',{ l x[})K^{X')T^-\-X',{ l x' i ) (66) 

with the grading ei = e 2 = 1. Here, we denote A' — X+^,fi' — fi+^. The reflecting matrix can also be 
interpreted as an operator matrix with higher spin. Explicitly, with the help of (0,p3), we have 

sin(2X' -n) ( A(X',c') B(X',c') 



K^(X') 



o i v: 



sin(2X') \ C{X',d) D{X',c') 



A [A) ~ \C+(X' -2) D+(\'-Z))> [ > 

where d = c+ ^. Note that the solution of the reflection equation can be changed by a gauge transforma- 
tion. In order to prove that the above defined nested transfer matrix is still a transfer matrix with higher 
spin reflecting matrix, we should prove that K^(X') and K^ (A') satisfy the reduced reflection equa- 
tion and its corresponding dual reflection equation. Indeed, it can be shown that the following reflection 
equation holds, 

f 12 (A' - ^)K[ 1] (A')r 2 i (A' + n')K^ (,/) = K™ (/i')f 12 (A' + M ')#f > (X')f 21 (A' - ft'). (68) 

With AfW = diag.(e 2lr] , 1), and the isomorphism (p2|), we find that K^- 1 ' satisfies the following relation 

f 12 (-A' + M ')^ (1) i h (A') stl M (1) r 1 ^i(27 7 -A'- M ')^ (1) 2 + (M') st2 M (1) 2" 1 

= ^ 1 ) 2 + (//)^Af( 1 ) 2 _1 f 12 (2r;-A'-//)^ 1) (A')Af( 1 )7 1 f 21 (-A' + / /). (69) 

By use of the cross-unitarity relation r^ 1 (2?7 — A)M{ f^i (A) M{ = sin(X)sin(2r] — A) • id., the above 
relation is just the dual reflection equation which we need. 

The row-to-row monodromy matrix T*- 1 ' (A', {/4}) (corresponding to the periodic boundary condition) 
and its inverse are defined as 

T^JX'A^Z-Z = ^(A / + Mi)^r^' + ^)a?S--r(A' + M' 1 )SI^ (70) 

t (1 C(-a', M})&::.£ = W - aO"" • • ■ M*' - /4)££MA' - /4)&v ( 71 ) 

We show that a problem to find the eigenvalue of the original transfer matrix £(A) reduces to a problem 
to find the eigenvalue of the nested transfer matrix fW (A). The nested transfer matrix is still a boundary 
case with higher spin reflecting matrix. 

12 



In order to ensure the assumed eigenvector is indeed the eigenvector of the transfer matrix, /ii, • • • , ji n 
should satisfy the following Bethe ansatz equations, 

sin(2u-i — Ti)sin{iLj — c + ?; — sri)sin(ui — c + 2r\ + sri) OM , 

V-^t: ; \ zr 1 - ; , ~ ~ , ' Sin 2N {jlj + rj) 

sin(2fij + rj)svn{jij — c + rj + srijsinyjjij — c — srj) 

a 

x []sm(/ij + fii)sin(fij - ji l - rj) = -sin 2N ' {pj)Ay> '(fij), j = 1, 2, • • • , n. (72) 

8=1 

Here we have used the notation AW(A) to denote the eigenvalue of the nested transfer matrix t^\\). 

4.2 Bethe ansatz for the six-vertex model with higher spin reflecting matrices 

We repeat almost the same procedure as that of the first level algebraic Bethe ansatz method. We only 
write down some results without the detailed calculations here. We have 

c «^ ^"(2^-77) > = c _ ax sin(2X)sin(X + C + V ~ sr,) > ^ ^ > ^ ,_ 

sin{2X') sin(2X + Tj)sin(X — c+ srj) 



n ^(2X'-v) [A{ X',c') + D(X\c') S : n ^r J |0 > f , (74) 

sin(2X') sin(2X — rj) 



i{2X+n) sin ( X + c + r] + sr))sin(X - c - rj + srj) > ^ > 

sin(X — c — r) — srj)sin(X — c + srj) 



(75) 



and 



A+(X> - 2)| > r = _ e ^ + ,) ""(A + c^ + ^) |0 ^ +|Q ^ (?6) 

2 sin{\ — c + rj + srj) 

r _i_, / V . 1,., 77 , sin(ri)e~ t2X ,, 

i2A sm(2A-T?)sm(A + c-?7-577)sm(A-c + ?7-gT?) |n ^ _ + 
s7n(2A)sin(A — c — srj)sin{X — c + 77 + S77) 

We write the nested double-row monodromy matrix as 

T (A ' {M) - (, C«(A) D«(A) ■) (?8) 

From the results obtained above, we know that this double-row monodromy matrix also satisfies the 
reflection equation (pq). Considering the transformation 

A^(X)=A^(X)- ^f^- V^(X), (79) 

siti(2A — rj) 

we can find the following commutation relations which are useful for the algebraic Bethe ansatz method, 

D«(A)CW(m) = ^-M + ^n(A + /i) c(1? (1) 

sin(X — fi)sin(X + ji — rj) 

- Sl f" >^]f-"\ C (D (A) P (D (,) + ""ff <( A+M \ C W (A )^(D (,), (80 ) 
s«n(A — ji)sin{2ji — rj) sinyX + ji — rj) 
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i( i) (A)c (i) (M) = «n(A- / x-^n(A + / x-2^) c(1)(/i)J[(1)(A) 

SWl(A — flJSlTlyX ~1~ fJ> — Tf) 

sin(X — /j,)sin(2X — 77) 
_ «n(2^n(2A - 2,) ri n(,)e-^) 
szn(A + /z — rj)sin{2X — r])sin(2^i — rj) 

C (1) (A)C (1) (/i) = C (1) (/^)C (1) (A). (82) 

We thus finally obtain the eigenvalues of the nested transfer matrix as 

a«<y) . -upa n W + w - rin fl { "" (A ;:,"' i 'lm:"";^ 5 " ,"' | 

i=i /=i l sm(A'-^'J ; )szn(A' + ^'J -77) J 

n 

-UtU 2 IJ[«n(A' + M'i - »?)a*n(A' - ^ - »?)] 

" / sm(A'- M f ) +7 ? ) a m(A' + / /< 1) ) | 
f=i I **»(*' - m'P W A ' + m'j (1) -»?)/' 

where /j ^ , • • • ,(j! m should satisfy the corresponding Bcthe ansatz equations. In waht follows, we give 
a summary of our main result. 

4.3 Result 

The eigenvalues of the transfer matrix for the generalized supersymmctric t — J model are given as follows: 

sin(2A — rf)sin{X — c + rj — srj)sin(X — c + 1r\ + sr\) 



A(A) = 



xsm 



sin(2X + r))sin(X — c + 77 + sr])sin(\ — c — sr/) 
2N/\ , .jj sin{\ + Hi)sin{\- Hi-rj) 



f (A + r?)n 



*■ sin(X + Hi + rj)sin(X — fii) 



1 
-sin 2N {X) J] -^-tt , ■ n , , , A (1) (A), (84) 



/jn sin(X + c + 77 + sr])sin(X — c— r/ + srj)sin(X + c — 77 + sr/) 

sin(X — c — r\ — sr\)sin{X — c + srj)sin(X — c + r\ + srj) 

x n^A + * + .)«n ( A - mo] ft j sm(A 7 ;(1) ~ (1 f m ! A + ill; ^ 

»=i (=1 I sin(X-n\ ')sin{X + Hi') 

sin(2X — rj)sin(X + c — 77 — sri)sin(X — c + 77 — sr])sin(X + c + 77 — st/) 
sm(2A + r/)sin(X — c — sr))sin(X — c + 77 + sr])sin(X — c + S77) 

sin(X — n\ \ + rj)sin(X + /j,\ + 77) 



x 



n[«m(A + w)a i„ ( A - w - v)] n — "7 T V ■ (85) 
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where n\ , • • • , jjb n and fj\ , • • • , n m should satisfy the Bethe ansatz equations 

sin{iA + c + r\ + srj)sin(iA — c — r\ + si])sin(fi- + c — r\ + S7})sin(ji,j — c — srf) 
sin((j,j — c — 7] — sr])sin(jjL; + c + r\ — srj)sin(fij + c — r\ — srj)sin(fi, — c + r\ — srf) 
JL sin(jij + Hi)sin(jij — iM — rj) Ji sin(ny' - /i ; (1) + j])sin((ij + /ij + r/) 
i=i sin((J,j + (Mi + i])sin(fif ) - m) l=1<jtj sin{^f ] - /i, (1) - ji)sin(^ 1] + ^ - rj) ' 

j = l,---,m, (86) 



and 



ain(jij + c — r] — sr])sin(X + c + 77 — s?i) sin 2N (/j,j + rj) -pr sin(nj — fi\ ')sin(jj,j + n\ ) 



n 



s™(ft - c + 277 + srj)sin{\ - c + arf) sin 2N {^) f = \ s i n ( Mj . _ ^W + ^) sm ( ft + ^W + r/ ) ' 

j = !,-■■, n. (87) 

5 Summary 

In this paper, we have studied the generalized supersymmetric t — J model with Kondo impurities in 
the boundaries. Using the higher spin L operator of XXZ Heisenberg chain and the general diagonal 
solution to the reflection equation for six vertex model, we find a higher spin reflecting matrix for the 
generalized supersymmetric t — J model. Applying the graded algebraic Bethe ansatz method, we obtain 
the eigenvalues of the transfer matrix for the t — J model with higher spin boundaries. 

It is interesting to solve this problem in other background gradings, for example, FBF or BFF. 
The higher spin reflecting matrix should be constructed from the BF or FB six vertex models. The 
analysis of ground state properties, low-lying excitations and thermodynamic Bethe ansatz is always 
worth performing. 

One can find that the SU q (2) higher spin reflecting matrix also satisfy the reflection equation of 
SU q (N) model. The eigenvalues of SU q (N) model with SU q (2) higher spin boundary impurities can be 
obtained by using the nested algebraic Bethe ansatz method. Actually, the SU q (2) higher spin boundary 
impurities could be embeded into SU q (M\N) spin chains with M >= 2 or N > 2. 

After we put our article to the cond-mat e-print archive, X.Y.Ge and H. Q.Zhou inform us that they 
solve the same problem independently |46j. 
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